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PROBING THE STATISTICS OF THE TEMPERATURE-DENSITY RELATION OF THE IGM 
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ABSTRACT 

Gravitational instability induces a simple correlation between the large and small scale fluctuations 
of the Ly-a flux spectrum. However, non-gravitational processes involved in structure formation and 
evolution will alter such a correlation. In this paper we explore how scatter in the temperature-density 
relation of the IGM reduces the gravitationally induced scale-scale correlation. By examining whether 
or not observations of the correlation are close to that predicted by pure gravity, this puts constraints 
on the scatter in the p — T relation and in turn on any physical process which would lead to scatter, 
e.g. strong fluctuations in the UV background or radiative transfer effects. By applying this method 
to high resolution Keck spectra of Q 1422+231 and HS 1946+7658, we find the predicted correlation 
signal induced by gravity, and the diminishing of this correlation signal at small scales. This suggests 
extra physics affects the small-scale structure of the forest, and we can constrain the scatter in the p — T 
relation to a conservative 20% upper limit. A crude model suggests, if there is any spatial correlation of 
temperature, the coherence length scale must be smaller than ~ 0.3/i _1 Mpc to be consistent with the 
Keck data. 

Subject headings: large-scale structure of the universe — methods: statistical — methods: numerical 
— intergalactic medium — cosmology: theory 



1. INTRODUCTION 

In the last few years our understanding of small-scale 
structure at high redshift has increased enormously through 
the study of the Ly-a forest, the series of absorption features 
in the spectra of distant QSOs. Based in large part upon 
insights gained from numerical simulations (see, e.g., Cen, 
Miralda-Escude, Ostriker, & Rauch 1994; Zhang, Anninos, 
,& Norman 1995; Hernquist, Katz, Weinberg, & Jordi 1996) 
we now know that the physics governing the Ly-a forest is 
relatively simple. At the relevant redshifts the gas making 
up the intergalactic medium (IGM) is in photoionization 
equilibrium, which results in a tight density-temperature 
relation for the absorbing material with the neutral hydro- 
gen density proportional to a power of the baryon density. 
Since pressure forces are sub-dominant, the neutral hydro- 
gen density closely traces the total matter density on the 
scales relevant to the forest (0.1 — 10/i" 1 Mpc). The struc- 
ture in QSO absorption thus traces, in a calculable way, 
slight fluctuations in the matter density of the universe back 
along the line of sight to the QSO, with most of the Ly-a 
forest arising from over-densities of a few times the mean 
density (see, e.g., Croft, Weinberg, Katz, & Hernquist 1998; 
Croft et al. 2002). 

While most attention has focused on the low-order statis- 
tics of the Ly-a forest, the statistical properties of the flux 
are significantly non-Gaussian. Viel et al. (2004) studied the 
bispectrum of a large sample of QSO Ly-a flux spectrum. 
While their results from observations and simulations and 
theoretical predictions reasonably agree with each other at 
large scales, they found the errors on the bispectrum are 
too large to discriminate between models with different Ly- 
a forest distribution based on the bispectrum method. 

However, while in general probing high-z physics with the 
bispectrum is complicated, in some special cases patterns of 
higher order correlations can be visible in the moments of 
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the flux. One particular pattern is that due to the action of 
gravity in hierarchical structure formation. The presence or 
absence of this pattern provides us with a powerful probe of 
'extra' physical effects, beyond gravity, in the formation of 
the Ly-a forest, as first emphasized by Zaldarriaga, Seljak 
& Hui (2001; hereafter ZSH01). These authors noted that 
it was possible to constrain fluctuations in the continuum 
by appealing to higher order moments of the flux. In this 
paper we wish to address another physical effect: scatter 
in the density-temperature relation induced for example by 
different thcrmodynamical histories for the gas, fluctuations 
in the UV background radiation field or radiative transfer 
effects. We shall see that it is possible to put tight con- 
straints on the scatter in the p—T relation by considering 
the higher order moments of the spectra. 

The outline of the paper is as follows. In §2 we briefly 
discuss the simulations involved in this work and how we 
extract "mock" Ly-a spectra. The statistical method intro- 
duced by ZSH01 is reviewed in §3. Our main results are 
presented in §4, and we discuss the application to high res- 
olution Keck spectra in section §5. In §6 we discuss several 
related issues. 

2. SIMULATION 

Since we are interested in the higher order correlations in 
the Ly-a flux which are induced by gravitational processes, 
we first need a model of gravitational clustering. For this 
purpose we use a particle mesh code to follow the evolution 
of dark matter clustering, including only gravity. Details 
of the PM simulation can be found in Meiksin, White, & 
Peacock (1999), White (1999) and Meiksin & White (2003). 
Given the particles' positions and velocities, the Ly-a forest 
spectrum can be constructed based on the so-called "fluc- 
tuating Gunn-Peterson approximation" , in which Ly-a ab- 
sorbing gas is assumed to trace the dark matter distribution, 
and be in a state of photoionization equilibrium with the 
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Fig. 1. — Statistics of the Ly-o flux spectrum based on a square filter 0.2 < k < 0.3. Solid lines show the no scatter case, dashed lines a = 0.1 
and dotted lines a = 0.2. Panel (a): power spectrum; panel (b): bispectrum; panel (c): trispectrum; and panel (d): cross-correlation coefficient 
C(k). In panel (d) C(k) is smoothed with a boxcar method for visual clarity. 



UV background radiation (see, e.g., Cen, Miralda-Escudc, 
Ostriker, & Rauch 1994; Zhang, Anninos, & Norman 1995; 
Hcrnquist, Katz, Weinberg, & Jordi 1996; Hui & Gncdin 
1997). Under these assumptions the gas temperature can 
be obtained through 



T = Tn 



(1) 



where p is the matter density and p the cosmic mean. The 
index 7 is expected to vary from near unity right after reion- 
ization, to <~ 1.5 at low-z (Hui & Gnedin 1997). We will set 
T = 2 x 10 4 K throughout, but in the next section we will 
introduce scatter around the values given by Eq. 1. 

The neutral hydrogen number density, based on our ion- 
ization equilibrium assumption, is nni oc (p/p) 2 T~ - 7 , with 
the constant of proportionality dependent on the baryon 
density, redshift, and hydrogen photoionization rate. We 
shall fix it using the observed value of the mean flux (see 
below). With the gas density, temperature and velocity 
data, the optical depth r, at a given velocity u is 



J nm(x) <y(x) dx 



dx 



(u - u ) 2 
b 2 



(2) 



where uo is the sum of Hubble expansion and peculiar ve- 
locity at position x, and b = ^j2ksT /m p is the Dopplcr 
parameter. Here m p is the proton mass, and a is the pho- 
toionization cross section. The flux at velocity u is then 
cxp(— t). 

For definiteness we concentrate on the z — 3 output of 
a 512 3 particle and 1024 3 mesh PM simulation in a box 
of side 60ft, -1 Mpc, which gives a simulation resolution of 



~ 6kms 1 at z ~ 3. The cosmology simulated was of the 
ACDM family, with fi mat = 0.3, h = 0.67 and <j s = 0.9. 
The density and velocity fields were Gaussian smoothed by 
lOOft^kpc to simulate the effects of gas pressure and 32 2 
spectra, each of 1024 pixels, were computed on a regular 
grid parallel to the box axes. The normalization constant 
was fixed to reproduced the mean flux, F = 0.684, observed 
by McDonald et al. (2000) at z = 3 (see also Appendix B 
of Meiksin & White 2004). 

3. STATISTICS 

To constrain scatter in the temperature density relation- 
ship on small scales we adopt the high-order statistical 
method proposed by ZSH01. We briefly describe this tech- 
nique here and refer the reader to their paper and Mandel- 
baum, McDonald, Seljak, & Cen (2003) for details. 

We define the relative fluctuation of the flux, Sf = 
F/F — 1, and its Fourier transform 



F(k) 



dxe lkx 5 F {x) 



(3) 



We bandpass filter 5f with a window function W{k): 
Ff{k) — F(k)W(k), taking the window function W(k) to 
be either a square window, 



W klM {k) = 



if k\ < k < ki\ 
otherwise 



(4) 



or a Gaussian window function. Our results are not sensi- 
tive to this choice, so in what follows we shall concentrate 
on the square window. 

The filtered field Ff(k) in Fourier space is then inverse 
transformed back into real space to define 



tkx 



Ff(k). 



(5) 
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Fig. 2. — Cross-correlation coefficient, C(k), for different filtering schemes. Symbols are the same as those in Figure 1. Panel (a): 
0.02 < k < 0.06; panel (b): 0.1 < k < 0.2; and panel (c): 0.2 < k < 0.3. Keck data is plotted as the solid circle. Error bars are 1<t and 
plotted every three points for visual clarity. 



and this is squared to produce h(x) = Sj I (x). Following 
the conventions in Mandelbaum, McDonald, Seljak, & Cen 
(2003), we define P as the power spectrum of Sp, T as the 
power spectrum of h, and B as the cross-spectrum between 
Sp and h. Since h depends quadratically on Sp, B and 
T can be related to the flux bispectrum and trispectrum, 
respectively. We also introduce the same cross-correlation 
coefficient C(k) to study the effect of large scale density 
fluctuation on small scale power, defining C = B/\JPT . 

4. SCATTER IN p - t 

We model the scatter in the temperature-density relation 
as log-normal in T at fixed p, with mean value determined 
from Eq. 1. We choose the width, <r, to mimic the results 
seen in hydrodynamic simulations. Specifically we choose 
a = ShiT = O.f and a — 0.2 to study the effect. For com- 
parison, by checking p — T scatter plot, we find a scatter 
of a = 0.1(0.2) is consistent with results from low (high) 
resolution hydrodynamic simulations in Bryan, Machacek, 
Anninos, & Norman (1999), respectively (see their Fig. 3). 
Figure la shows the power spectrum P of our simulation 
under the assumptions of no scatter, a = 0.1 or a = 0.2. 
At k > 0.1 skm -1 (scales < 10 /i~ 1 Mpc), scatter in the 
temperature-density relation adds power. 

We show in Fig. 1 the power spectrum, bispectrum, 
trispectrum and cross-correlation coefficient from the simu- 
lations with a top-hat filter with modes between 0.2 < k < 
0.3 skm^ 1 . We chose that range to maximize the effect of 
p — T scatter on C(k). As pointed out by ZSH01, a posi- 
tive correlation between the large scale density fluctuations 
and the power on small scales is induced by gravitational 
instability. Since, within the "fluctuating Gunn-Peterson 
approximation" , the fluctuations in the dark matter distri- 
bution are directly traced by fluctuations in the Ly-a forest, 
the positive correlation induced by gravitational instability 
is manifest in the Ly-a flux. In the Ly-a forest the cor- 
relation is negative since the flux decreases with increasing 
density. 

We see that our model for scatter in the p — T relation 
does not change the large scale fluctuations of the Ly-a flux 
spectrum; however, by locally varying the optical depth t 
through Eq. 2, it effectively adds noise at small scales. Fig- 
ure Id shows C(k) at low-A: becomes steadily less negative as 
the scatter is increased. Note also that the bispectrum, tri- 
spectrum and cross-correlation coefficient go quickly to zero 
for k > 0.1 skm -1 . This is because after filtering the high- 



est mode in the squared field, h(k), is k 2 — k\, or 0.1 skm -1 
in this case. 

Finally, we have investigated the requirements on signal- 
to-noise and resolution of spectra which would be suitable 
for measuring C(k) on the appropriate scales. By adding 
Gaussian random pixel noise to our mock spectra we found 
that a S/N approaching 50 - 100 is required in order not 
to wash out the small scale signal. We also notice that 
the level of scatter is closely related to the simulation res- 
olution. The flux spectrum is obtained by averaging out 
fluctuations in simulations cells over the thermal broaden- 
ing width. The smaller the cells are, the more cells will 
be averaging out, and the smaller the fluctuation will be. 
In our simulation, the typical temperature of the IGM is 
~ 2x 10 4 K, which corresponds to a size of <~ 18kms _1 . 
So the flux spectrum in each pixel is roughly averaged over 
three neighboring cells. To test this effect, we check simu- 
lations with two other resolutions that are close to the ones 
in hydrodynamic simulations by Bryan, Machacek, Anni- 
nos, & Norman (1999), which are two and four times higher 
than the one used in the paper. We found that there are no 
significant differences among these cases. 

5. AN APPLICATION TO QUASAR DATA 

In this section, we study the cross-correlation coefficient 
C{k) with two, high resolution, high S/N quasar spec- 
tra, Q 1422+231 at z = 3.62 (Rauch, Sargent, Barlow, 
& Carswell 2001) and HS 1946+7658 at z = 3.05 (Kirk- 
man & Tytler 1997). These two quasars were observed with 
Keck High Resolution Echelle Spectrometer (HIRES), with 
a signal-to-noise ratio (S/N) of <~ 140 and ~ 50 — 100, 
respectively. We split the two spectra into a total of 13 
separate pieces between Ly-a and Ly-/3 emission to avoid 
contamination from Ly-/3 absorption. Each piece has a co- 
moving length of 60/j" 1 Mpc to match our PM simulation. 
We calculate C(k) among these 13 samples following the 
method described in section §3, and compare these results 
against our PM simulation. We made no attempt to correct 
for metal lines. In order to account for evolution in the mean 
flux, we adjusted each spectral segment to (F) = 0.684 by 
scaling the implied optical depths as in the simulations. 

In Figure 2 we show C(k) from both Keck data and our 
PM simulation. To study the filter effect, we select three 
filters: (a): 0.02 < k < 0.06 skm" 1 ; (b): 0.1 < k < 
0.2 skm" 1 ; and (c): 0.2 < k < 0.3 skm" 1 . With each fil- 
ter we plot the no scatter case (solid lines), scattering with 
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a = 0.1 (dotted lines) and scattering with a = 0.2 (dashed 
lines) . Keck data are labeled as solid circles, with error bars 
computed from the simulations (the sample size in the data 
is too small to give constraints on variances). 

Clearly, with large fluctuations, negative signals present 
in cases with filter (a) & (b), which agrees with simulation 
quite well. This is consistent with what ZSH01 found in 
Q 1422+231 data, and confirms the correlation between the 
large and small scale fluctuations of the Ly-a flux spectrum 
induced by gravitational instability. Also in case (a) & (b) , 
there is almost no difference between no-scatter and scatter 
because their power spectra (see Figure la) show no dif- 
ference between 0.01-0.2 skm^ 1 . At smaller scales (larger 
k modes) the scattering increases the power spectrum and 
starts to differentiate between the scatter and no-scatter 
cases. We also find, with filter (c), the C(k) from Keck 
data is essentially zero, consistent with 20% scatter in p — T 
relation case. Though other effects (such as random noise) 
may also drive C(k) to zero, increasing the p — T scatter 
above 20% drives C(k) positive, so we can conservatively 
set an upper limit of 20% on the the scatter in the p — T 
relation. 

In the above we adopted an independent random 
Gaussian-distributed T for each pixel. In reality, however, 
pixels tend to correlate with neighbor pixels: pixels around 
a high T pixel appear to also have a high temperature and 
vice versa. Such a correlation among pixels tends to pro- 
duce a "correlated" noise from p—T scatter and enhance the 
signals in cross-correlation coefficient C(k): the larger the 
correlated pixel length scale, the stronger the C(k). So with 
filter (c), we can ask: what is the largest length scale that 
pixels can be correlated with each other, without showing a 
strong signal in C(k) that conflicts with the Keck data? We 
can parameterize this by adopting a correlation parameter 
r: Tj = rTi_i + (1 — r)T now , while Tj is the temperature 
of pixel i, T ncw is the temperature obtained from an inde- 
pendent random draw, and r is a number between and 
1: r = means pixels are independent of each other, while 
r = 1 mean the maximum correlation. We found, for the 
20% scatter case, r must be smaller than <~ 0.3 - 0.5 not to 
be ruled out by the Keck data. This roughly corresponds 
to a distance of <~ 5 pixels if we define two pixels correlated 
at the 1% level as "correlated". This means, with a 20% 
scatter, the correlated scale length must be smaller than 
- 0.3ft,- 1 Mpc, or - SOkms" 1 . 

The observational evidence for the gravitationally in- 
duced correlation, and for small-scale decorrelation, cur- 
rently rests on two high resolution spectra. It would obvi- 
ously be of some interest to improve the analysis presented 



here and to increase the sample of QSO spectra for which 
the analysis has been performed. 

6. DISCUSSION 

ZSH01 showed that gravitational instability, and the as- 
sumption of the "fluctuating Gunn-Peterson approxima- 
tion", predict a specific correlation between large- and 
small-scale power in the Ly-a forest. They also showed, 
based on high resolution Keck data, that observations were 
consistent with the correlation predicted by gravity. Here 
we point out that this observation allows one to limit the 
deviations from the "fluctuating Gunn-Peterson approxima- 
tion" and in particular to constrain the fluctuations in the 
density-temperature relation of the IGM. Such fluctuations 
are expected to arise from varying thermodynamical his- 
tories for the gas and from fluctuations in the ultra-violet 
background radiation field or radiative transfer effects. We 
find that these fluctuations must be < 20% in T at fixed p 
in order to be consistent with high resolution Keck data. 

Our method provides a sensitive way to probe for scatter 
in the temperature density relation that is expected from 
the complex physics of the IGM and seen in cosmological 
hydrodynamic simulations. Persistence of gravitationally 
induced higher order correlations in the spectra requires 
that the scatter be no more than twice as large as already 
seen in hydrodynamic simulations of structure formation. 
This limits the role of "additional" sources of fluctuations 
such as an inhomogeneous UV background field or radiative 
transfer effects. To further strengthen the case requires this 
correlation to be computed from more high S/N, high reso- 
lution quasar spectra such as could be obtained with Keck, 
Magellan, or similar telescopes. We hope that when more 
high quality spectra along different lines of sight become 
available, this method will be a powerful tool to study any 
"extra" physical process, beyond gravity, in the IGM. 
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